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Analytical models are presented for calculating mass flow and 
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orifice or a converging-diverging nozzle. The solutions for both 
frozen and shifting equilibrium flow assumptions give comparable 
results, yet overcome difficulties encountered with previously used 
models. These solutions, purely analytical, are applicable for any 
liquid-vapor fluid mixture and require only nozzle inlet conditions for 
making calculations. 
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DEFINITION OF SYMBOLS 


Symbol 

A 

B 

c F 

C c 

F 

G 

g 

H 

I 

J 

P 

Q 

R 

S 

T 

U 


Definition 

Passage cross sectional area (ft 2 ) 
Ratio of gas to liquid velocity, Ug/Uj 
Nozzle thrust coefficient 
Contraction coefficient, A c /A^- 
Thrust (lbf) 


Mass flow rate per unit area, W/A ^~~^ z ) 

( ft lb m \ 

7 1 ; ) 

sec^ lbf / 


Enthalpy (Btu/lbm) 

„ ... . . /lbf sec\ 

Specific impulse I — — — — I 


Mechanical equivalent of heat (778) 



Static pressure 



Fluid mixture quality, Wg/W m 


Gas constant 


ft lbf \ 
lb m 0 R / 


/ Btu \ 

Entr °P' r 

Absolute temperature (°R) 
Velocity (ft /sec) 


v 



Symbol 

W 

Y 

P 

£ 

4 > 


DEFINITION OF SYMBOLS (CONT. ) 

Definition 


Mass flow rate 



Specific heat ratio, Cp/C v , of vapor phase at inlet plane 


Fluid density 



Nozzle expansion ratio, A x /A^ 

Two-phase critical flow parameter, equation (16) 


Subscripts 

a 

c 

e 

g 

i 

l 

m 

t 

x 

Ho 


Ambient 

Vena contracta plane 
Nozzle exit plane 

Designates mass of fluid that is gas at the inlet plane 
Nozzle inlet plane 

Designates mass of fluid that is liquid at the inlet plane 
Total two-phase fluid mixture 
Nozzle throat plane 
Station along nozzle axis 

Designates calculations were based on homogeneous 
shifting equilibrium conditions 


vi 



TECHNICAL MEMORANDUM X-53269 


CALCULATIONS OF MASS-FLOW AND THRUST 
PRODUCED FOR A TWO-PHASE FLUID 
MIXTURE PASSING THROUGH A CHOKED NOZZLE 


SUMMARY 


Analytical models are presented for calculating mass flow and 
thrust produced for a two-phase fluid mixture passing through a choked 
orifice or a converging-diverging nozzle. The solutions for both 
frozen and shifting equilibrium flow assumptions give comparable 
results, yet overcome difficulties encountered with previously used 
models. These solutions purely analytical, are applicable for any 
liquid-vapor fluid mixture and require only nozzle inlet conditions for 
making calculations. 


INTRODUCTION 


Cryogenic propellant tanks coasting in space must be vented 
periodically to prevent overpressurization. Since the vehicle is under 
a zero or near-zero gravity condition, it is possible for a two-phase 
fluid mixture to enter the vent system and be discharged overboard. 
The variations of flow and thrust produced when such a mixture passes 
through the discharge ports must be considered in the vent system 
design. 


From a literature survey made by R. V. Smith (Ref. 4), it was 
found that no satisfactory solution exists for calculating choked two- 
phase flow rates for the full range of fluid qualities. FIGURE 1 shows 
the available experimental data and the results of various analytical 
calculations surveyed by Smith. It is evident from this figure that 
none of the reported analytical solutions are applicable for the entire 
fluid quality range. Other difficulties, such as the use of empirical 
constants, requirements of known throat conditions, and difficulties 
of incorporating thrust calculations, make the use of the reported 
analytical solutions undesirable. The methods for calculating two- 
phase flow described herein overcome the above difficulties and give 
comparable results with previously used models. 


It can be seen from FIG 1 that a considerable variation in 
experimental mass flow occurs for a given fluid mixture quality. This 
flow variation is greatest for low quality fluid conditions, therefore, 
it is thought to be caused by the rate at which the fluid properties shift 
toward equilibrium conditions during the expansion process. In an 
attempt to account for this flow variation, two approaches for calcu- 
lating mass flow are taken in this analysis. The first flow model 
assumes that the fluid properties remain under frozen equilibrium 
conditions during the expansion process while the second model assumes 
that the fluid properties immediately shift to equilibrium conditions. 

The flow models described herein should predict the upper and lower 
mass flow limits and the thrust produced when a two-phase fluid mixture 
passes through a choked orifice or a conver ging-diver ging nozzle. 


THEORY 


General Flow Equations 


The following basic equations are used in calculating the two- 
phase mass flow and nozzle thrust for both the frozen-equilibrium and 
separate -phase shifting equilibrium flow models and are given to define 
clearly the symbols used in this report. The quality, Q^, of the two- 
phase fluid entering the nozzle is the ratio of vapor to total mass flow 
and is defined by 


Q iS 


w 


SL 


w ii+ w gi 


(1) 


The flow passage cross sectional area, A x , at a given plane 
location, x, along the nozzle axis is equal to the sum of the cross 
sectional areas occupied by the liquid and vapor phases. 

^x = ^gx + -Aj x (2) 

From the continuity equation, the mass flow per unit area for 
each phase is 


Gi 
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By combining equations (1), (2), (3) and (4) and simplifying, 
the total mass flow per unit area is obtained. 


J mx 


1 - Q 


Gi: 


i 


Qj 


G 


g x 


(5) 


The total thrust produced by a nozzle with the exit plane at any 
station, x, equals the summation of the thrust produced by the separate 
phases. 


F mx = x + F gx (6) 

The thrust produced by each phase is calculated from the 
momentum equation, 

F lx= — * g W< * + (Pi x - Pa)Aj x (7) 

F gx = + <P gX - Pa)A gx (8) 

E> 


Since it is assumed that the static pressures of the separate 
phases along the nozzle wall are equal (Pj x = Pgx = Px)> equations (1), 
(2), (6), (7) and (8) can be combined and simplified to give the overall 
nozzle thrust coefficient for the two-phase fluid mixture. 
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where 
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Knowing the thrust coefficient, one may calculate the specific 
impulse for the two-phase fluid. 
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( 11 ) 


Cp Pj 

T r mx 1 

-*-mx ” p; 

Calculations for Frozen Equilibrium Model 

The following assumptions describe the nozzle or orifice flow 
system and are implied when referring to "Frozen Equilibrium 
Conditions. 11 


a. The system is frictionless with no interaction or energy 
exchange between the liquid and vapor phases. 


b. Metastable, isentropic, one-dimensional expansion occurs, 
i. e. , no evaporation or condensation occurs during the expansion 
process and the vapor phase acts as a perfect gas. 


c. The static pressures of the separate phases are equal at 
any plane along the nozzle axis. 

d. The fluid velocity at the nozzle inlet plane equals zero. 


The velocity of the liquid and vapor phases at any station along 
the nozzle axis may be calculated from the general energy and perfect 
gas equations respectively. 


UixV 2 g< P i ' P x)/Pii 



V - l 

(Px/Pi) ^ 


( 12 ) 


(13) 


By combining equations (12) and (13) with the continuity equation, 
the mass flow per unit area at any nozzle station can be calculated* 
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The nozzle throat or maximum flow conditions can now be 
determined by combining equations (14) and (15) with (5), differentiating 
G mx with respect to P x /Pi and equating to zero. 


dG mx 

d(p x / Pi ) 


(■-ft) 

yfv+i' 

\ V > 

it*) 

Si! 

2 - y 

1 * 

-(ft)'] 

(-7 

0.5 

h) 



S) 

Y+1-, 

Y 

1. 5 

(16) 


where 



Using equation (16), the critical two-phase flow parameter, i|j, 
was calculated and plotted on FIG 2 versus critical pressure ratio, 

Pt/Pi> for various vapor phase specific heat ratios, y. 

Knowing fluid properties at the nozzle inlet plane, one can 
readily calculate 4* and read Pt /Pi from FIG 2. After obtaining P^/P^, 
the critical mass flow, G m -t-, can be calculated using equations (5), (14) 
and (15). Substituting fluid inlet properties and G m t into the necessary 
equations (5 through 15), one can now calculate the nozzle thrust coeffi- 
cient, expansion ratio, and specific impulse at selected pressure 
intervals in the expansion process. In some problems where nozzle 
thrust must be calculated at various nozzle expansion ratios for different 
inlet qualities, it may be desirable to use a graphical solution similar to 
the one employed in the separate -phase and homogeneous models for 
calculating nozzle performance. 

Calculations for Separate -Phase Shifting Equilibrium Model 

The following assumptions describe the nozzle flow system and 
are implied when referring to "Separate-Phase Shifting Equilibrium 
Conditions. " 


a. The system is frictionless with no interaction or energy 
exchange between the inlet liquid and vapor masses. 
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b. Isentropic one -dimensional expansion occurs with fluid 
properties of the individual phase shifting to equilibrium conditions at 
all points during the expansion process. 

c. The portion of liquid that vaporizes or gas that condenses 
during the expansion process moves at the same velocity as its respec- 
tive initial phase. 

d. The static pressures of the separate phases are equal at 
any plane along the nozzle axis. 

e. The fluid velocity at the nozzle inlet plane equals zero. 

Since there are no simple equations for determining transport 
properties for a saturated fluid expanding along an isentrope, it is 
necessary to employ a graphical solution to simplify the critical flow 
and nozzle thrust calculations. The velocity at given points in the 
expansion process can be determined for both the liquid and vapor 
phases from the general energy equation. 


U ix =^gJ(H li '-H ix ) 

(17) 

Ug X =V^gJ(Hgi - Hg X ) 

(18) 


The mass flow rate per unit area for the separate phases can 
be determined from equations (17), (18) and the continuity equations 
(equations (3) and (4)). The values of enthalpy, H, and fluid density, 
p, for the separate phases can be obtained by starting from a known 
inlet condition on a pressure, P, versus entropy, S, chart and 
reading the required values at selected pressure intervals along an 
isentrope. The velocity and mass flow per unit area for both the liquid 
and gas phases can be calculated for the selected pressure intervals 
(FIGS 3 and 4). Using equations (3), (4), (5), (17) and (18), one can 
now determine the two-phase mass flow at selected pressure intervals 
* for a given inlet fluid quality. By plotting the calculated two-phase 
flow per unit area versus pressure (similar to FIG 4), the maximum 
mass flow for the mixture, G m t, and throat static pressure, Pt> are 
obtained. 

Knowing G m t and calculating Gmx from equation (10) one can 
select the static pressure, P x , from the mass flow plot (FIG 4) for 
desired expansion ratios. Using G m t> Px> Qi and selecting the phase 
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velocities from a plot similar to FIG 3, one can calculate the nozzle 
thrust coefficient and fluid specific impulse from equations (9) and (11), 
respectively. 


DISCUSSION OF RESULTS 


The greatest difficulty in making two-phase flow calculations is 
determining the ratio of liquid to gas velocities. In the homogeneous 
shifting equilibrium model (Appendix A), it is assumed that both the 
liquid and vapor travel at the same velocity. This assumption of equal 
velocity simplifies the solution, but the results do not agree well with 
experimental data; therefore, a velocity variation between the liquid 
and vapor appears to exist. To determine the separate -phase veloci- 
ties in both the frozen equilibrium and separate -phase shifting 
equilibrium models, it was assumed that the initial liquid and gaseous 
masses act independently and are accelerated by their own available 
energy. This assumption permits a purely analytical method for 
determining phase velocities for both the converging and diverging 
nozzle sections. 

Other assumptions used in describing the flow systems, such 
as the frictionless isentropic expansion, have produced results which 
are in reasonably good agreement with experimental data when used 
in making single phase critical flow calculations. Also, the one- 
dimensional expansion assumption is valid for properly designed 
converging -diverging nozzles and for well-rounded orifices. For a 
sharp-edged orifice, this assumption becomes less valid, and a contrac- 
tion coefficient, C c , should be used for determining the effective critical 
flow area. If a contraction coefficient is employed in determining 
critical mass flow, care should be taken to insure proper use of this 
coefficient in the thrust equations. 

The assumption of zero fluid velocity at the nozzle inlet plane 
was made to simplify the flow equations. In most instances the effect 
of inlet velocities on nozzle performance is negligible, however, a 
correction for inlet velocities other than zero may be made by using 
total rather than static fluid properties at the nozzle inlet plane. 

Since the homogeneous flow model assumptions permit the most 
straightforward solution to two-phase flow problems, it was used by 
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Smith (Ref. 4) and in this report as a basis for comparing the results 
of experimental and other analytical solutions. Critical mass flow 
calculations for various quality steam-water mixtures were made using 
the presented models, and the results are compared with experimental 
data on FIG 5. FIGURE 6 is a similar plot showing the calculated 
results for a two-phase hydrogen mixture. 

It may be noted on FIG 5 that the separate -phase shifting 
equilibrium model gives better agreement for the lower limit of experi- 
mental data than does the homogeneous model. This indicates that, 
when shifting equilibrium conditions occur, the mass flow should be 
calculated using the separate -phase rather than the previously used 
homogeneous model. For metastable flow conditions, Smith recommends 
using the vapor choking model for calculating critical mass flow. The 
results obtained from this model agree with results from the frozen 
equilibrium model for high quality-steam-water mixtures (FIGS 1 and 
5). However, for low quality mixtures and for fluids other than 
steam -water, the vapor choking and frozen equilibrium models do not 
agree. The ratio of mass flow calculated from the vapor choking model 
to that calculated from the homogeneous model is essentially the same 
for fluids such as water, hydrogen, oxygen, and nitrogen (Ref. 4); the 
results using the frozen equilibrium model vary with different fluids 
(compare FIGS 5 and 6). 

Until additional experimental critical two-phase flow data 
are available for fluids other than steam-water and for various nozzle 
configurations, it will be difficult to predict which theoretical model 
gives better results. Recent experimental data from the National 
Bureau of Standards (Ref. 1), indicate that saturated liquid hydrogen 
or liquid nitrogen expanding through a sharp-edged orifice follows the 
frozen equilibrium flow assumptions. However, other data (Ref. 2) 
obtained by expanding saturated gaseous hydrogen through a converging- 
diverging nozzle, indicate that shifting equilibrium conditions prevail. 
These experimental results show that either shifting or frozen equilib- 
rium flow may occur, depending on the fluid properties and nozzle 
configuration. 


8 



CONCLUSIONS 


Theoretical models presented give purely analytical solutions 
to two-phase critical flow problems. These solutions offer the 
following advantages over previously used methods. 

a. The solutions permit calculations to be made knowing only 
nozzle inlet conditions and fluid properties. 

b. The models allow calculations of not only critical mass 
flow but nozzle performance and fluid conditions at any point in the 
expansion process as well. 

c. The solutions are applicable for the entire range of fluid 
qualities. 


d. Although this report deals with two phases of a single 
fluid (i. e. , steam-water, GH 2 -LH 2 , etc.), the presented solutions 
may be equally applied to unlike fluids (air-water, GHe-LH 2 , etc. ) 
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APPENDIX 


HOMOGENEOUS SHIFTING EQUILIBRIUM MODEL 

Since the homogeneous shifting equilibrium assumption permits 
the simplest solution to two-phase flow problems, it is generally used 
as a basis for comparing results of experimental and other analytical 
methods. A number of analytical approaches have been proposed for 
determining the critical flow rate assuming homogeneous conditions 
(Ref. 4). Essentially the same results are obtained when using any of 
the proposed solutions, therefore, the best approach depends upon the 
available information. 

If calculations are to be made using the separate-phase shifting 
equilibrium model, it is recommended that the method described in 
this appendix be used for determining the homogeneous flow since the 
required data are the same. The following assumptions are implied 
when referring to "Homogeneous Shifting Equilibrium Conditions. " 

a. Isentropic one -dimensional expansion occurs, and the 
system is frictionless. 

b. The fluid mixture shifts to equilibrium conditions at all 
planes along the nozzle axis (free exchange of energy). 

c. The liquid and vapor phases travel at the same velocity. 

d. The fluid velocity at the nozzle inlet plane equals zero. 
Mass Flow Equations 


Knowing the inlet conditions, one can calculate the velocity of 
the mixture, U mx , at various pressure intervals along an isentrope. 


^mx = >/2gJ(H m i - H mx ) 


(A- 1) 


If the enthalpy of the liquid, H{ x , and the vapor, Hg X , phases 
are known, the enthalpy of the mixture may be calculated. 

W„ X H PX + Wi x H» x 

H mx =- J S T7T — = QiH gx + (1 - Qi)H ix (A-2) 


mx 
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Combining equations (A-l) and (A -2) and simplifying, one can 
obtain the mixture velocity from the individual phase velocities used in 
the separate -phase shifting equilibrium method. 

U mx = u ix^ 2 Q i + d " QiJ (A-3) 


where 


£5 

Uix 


Using the above equation and velocities calculated by the separate - 
phase shifting equilibrium model (FIG 3), one can calculate values of 
mixture velocity at various pressure intervals and mixture qualities. 
FIGURE A-l is a typical plot showing mixture velocities along the nozzle 
axis. The mass flow rate of the two-phase mixture is obtained from the 
continuity equation. 


"mx 


= U 


mx 


Pmx 


(A-4) 


By solving for the mixture density, p mx , i- n terms of the 
separate -phase velocities and mass flow rates, the equation for the 
mixture mass flow rate becomes 


Gmx 


Vb 2 Qj + (1 - Qj) 
1 - Qi I Qi B 
G i x G gx 


(A-5) 


A plot of mass flow rate versus pressure similar to FIG A-2 can 
be made by using equation (A-5). From these curves the critical mass 
flow rate of the mixture, G m t> can be obtained. 


Thrust Equations 


The thrust produced when a homogeneous two-phase mixture 
passes through a converging-diverging nozzle may be calculated from 
the following equations. 

W v-yi f U ■KY'I v 

F mx = + (Px - Pa)A x (A-6) 
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By modifying equation (A-6) the thrust coefficient may be 
calculated. 


_ _ Gmt U mx (Px " F a) e x 

F mx " Pi At ~ gPi + Pi 

where 

_ Ax _ G m t 
£X ~ = G^ * 


Using equation (A-7) and selecting required values from figures 
similar to (A-l) and (A-2), one can determine the thrust coefficient for 
various mixture qualities and nozzle expansion ratios. The specific 
impulse may now be determined. 


t mx 


C F Pi 
x mx 1 

^mt 


(A -8) 
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FIGURE A-l FLUID VELOCITY ALONG NOZZLE AXIS 
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